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Resonance vibrations of impact oscillator with biharmonic excitation 


Vladimir Burd 

Department of Mathematics, Yaroslavl State University, Russia 

Abstract. We consider a damped impact oscillator subject to the action of a biharmonic 
force. The conditions for the existence and stability of almost periodic resonance solutions 
are investigated. 


1 Introduction 

Representations of vibro-impact processes in terms of integrals of motion for vibro-impact 
system with one degree of freedom is determined by the physics of vibro-impact system. 
The corresponding variables (impulse - phase) are very natural for the vibrating processes 
accompanied by impacts. The impulse can be interpreted as the force characteristics of 
impact. The phase is the instant of collision. 

In case of an instantaneous impact the representations of periodic vibro-impact processes 
can be written via so-called periodic Green’s functions. Periodic Green’s functions are 
determined only by the linear part of the system. They appear below or in the form of 
Fourier series or hnite relations on intervals of periodicity. The impact force is represented 
by the Dirac 5-function and the values depend on the impulse and phase of the impact. The 
periodic Green’s functions are calculated using well known methods [2]. Similar approach 
was proposed by V.I. Babitsky and M.Z. Kolovsky [3]. 

We want to give a correct description of resonant modes in vibro-impact system with 
one degree of freedom. 

We chose to study a model with elastic impact and viscous damping. Such a model is 
discussed in [2], [14]. 

A model that takes into account the energy loss at impact is proposed in [2]. It is not 
clear how to explore this model on a mathematical level of rigor. 

Our study is based on a combination of the method of averaging on inhnite interval 
[6] and the methods of singular perturbation theory [10,13]. 

We will study the motion of conservative vibro-impact systems with one degree of 
freedom and subject to small biharmonic perturbations. This problem is reduced to study 
two-dimensional systems with fast rotating phase and slow varying coefficients. 

Periodic perturbations of smooth two-dimensional systems with fast rotating phase 
and slow varying coefficients were studied in [9]. The conditions of closeness of solutions of 
exact and averaged equations on a hnite asymptotically large time interval are established. 
Resonance almost periodic oscillations in such systems are investigated in [5] and [6]. 

Recently, much attention has been focused on studies of dynamics of nonlinear systems 
perturbed by a biharmonic external force with different frequencies (see, for example [1, 
4, 7, 8]). 
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2 Conservative impact oscillator 

Here we follow the method from the book [2], Consider a linear oscillator 

X + = 0 . 

At the point x = A, we arrange an immovable limiter and assume that once the 
coordinate x reaches the value A, an instant elastic impact occurs in the system so that 
if a; = A at the time instant ta, then the relation 

X(ta - 0) = -X(ta + 0) (1) 

holds. 

If A > 0 and the energy level in the linear system is insufficient to attain the level x = 
A, then linear oscillations with the frequency hi take place. In the presence of collisions, 
the oscillation frequency satishes the inequality ca > hi and rises as the energy rises but 
does not exceed the value 2VL. Hence 

VL<uj <2VL, a > 0. (2) 

If A < 0, the oscillation frequency u obeys the inequality 

212 < a; < oo, A < 0. (3) 

At A = 0 the image point passes any phase trajectory for the same time with the 
doubled velocity 2VL so that 

oj = 2VL, A = 0. (4) 

Condition (1) suggests that variation of the impulse <ho in the neighborhood of the 
impact instant ta takes the form 

J = x_ — x+ = 2x_, i;_ > 0, 

where x^ = x(ta =F 0). 

The resulting force becomes localized at t = ta- Hence 

^o\t=t^ = JS{t - ta) (5) 

and 

icK+0 

J <hodt = J. 

ta — O 

Impacts occur periodically when ta = to + aT, where a is an integer, and T is the 
period between impacts calculated by the equality T = 27ra;“^ and (2)-(4). Thus, for 
—cxD < t < oo, we obtain a T-periodic continuation of (5) 

$0 = JSrit — to). 
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where is the T-periodic ^-function. 

Solution of the equation 

x + + <ho(a;, i) = 0 (6) 

is understood as the T-periodic function x{t) such that its substitution into this equation 
transforms it into a correct equality (from the viewpoint of the theory of distributions) of 
the form 

X -I- D?x -f JSrit — to) = 0, 
where to is an arbitrary constant, and for all a = 0, ± 1 ,... 

x{to + aT) = A, J = 2x-{to + aT). 

At the same time, the restrictions 


x{t) < A, i;_ > 0 


are fulhlled, and the periods of oscillations, depending on the sign of A, £t the frequency 
ranges of (2)-(4). 

To describe the solution analytically, we assume to = 0. In this case, for 0 < t < Tq, 
the solution of equation ( 6 ) has the form 


X{t) = -jK[uJo{J){t - to),UJo{J)], 

^0 I ^0 cos kujQt T(/ 1 \ - 

27r02 TT 2^k=l 02_fc2^2 , - 


l^(t, Uo) 20 sin(Oro/2) 

—2r2Atan^^, J > 0, 


and the third relation here determines the smooth dependence u}o{J) at A 7 ^ 0 , whereas 
cuo = 212 at A = 0. 

Geometric conditions of an impact result in frequency intervals (2)-(4). Note that 
when A = 0 , the solution x{t) for 0 < t < tt/Q takes the form 


x{t) = -—smQt, 

where J is a frequency-independent arbitrary constant. 


3 Perturbed impact oscillator 

Now consider a perturbed impact oscillator 

X+ n‘^x+ ^o{x,x) = e[f{t,T)-'fx], (7) 

where e > 0 is a small parameter, r = et is a slow time, 7 > 0 is a constant. The function 
$ 0 ( 2 ^) i) describes the force of impact. Representation of this function is given in the 
previous section. 

We consider the perturbations of two types. First a biharmonic perturbation is the 
sum of two small periodic forces with close frequencies. The corresponding perturbation 
has the form 

fit, t) = ai sin ut + 02 sin(z/t -|- Tr), 
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where ai, 02 , T are real positive numbers. The perturbation function /(t,r) can be 
represented as 

/(t,r) = E(r)sin(z/t+ /3(r)), t = et, ( 8 ) 

where 

E{t) = A/af'T^aaa^"cosTV+^, cos /3{t) = 
sin/3(r) = 

The function ( 8 ) is periodic in t with the period 21: jv and periodic in r with the period 
27r/r. The function E{t) is strictly positive if Oi 7 ^ 02 , which will be assumed. Consequently 
the function f{t,et) is almost periodic function (see [ 6 ]) with two basic frequencies. 

Second perturbation is a biharmonic force with very different frequencies. This force 
has the form 

/(f, r) = A sin(z/t + 6 *) + i? sin Tr, r = et. (9) 

where A, B, u, T, 6 are real numbers. 

We assume ip = ujQ{J)t and transform equation (7) into the system in terms of the 
variables J,ip (impulse-phase), by making a change 


where 


k[iP,ujo{J)] 


X = -Jk[iP,ujo{J)], 

X = -JUo{J)K^[lp,Uo{J)], 


K{lp,j) 


UJo{J) ^ 


1 1 coskip 

271^“^ IT — k'^ 

'J k=l ^ 


( 10 ) 


hlo — 


n 

Uo{J) 


It follows from the theory of Fourier series (see, for example, [12, chapter 4]) that the 
function k^[iP,uq{J)] has a hnite discontinuities at the points ip = 2 / 7 r, where I is an 
integer, and continuous at all other points. Moreover the function K^[ip,uo{J)] has a 
derative at interior points of the intervals [27r/,27r(/ -|- 1)]. Fourier series of function 
K,^[ip,uo{J)] has a form 

1 k sin kip 
-MJpg) -2 ^q 2 f.2- 

k=l ^ 

Therefore change (10) is not smooth at '^ = 217t. Thus, in new variables, impacts occur 
when Ip = 21tt. Making the substitution (10), we arrive at the system 

^ = -Aeuo{J)[f{t,T) +-iJuQ{J)K^{ip,J)]K^{ip,J), 

f = Uo{J) - Aeuo{J)J-\f\kT) + Juo{J)k^{iP, J)](-J«(^, J))j. ^ ’ 

A detailed derivation of the system (11) is contained in the book [2, chapter 2]. This is a 
system with a fast rotating phase, where the right-hand sides are periodic in ip and have 
hnite discontinuities at the points ip = 21tt. The dependencies Ci;o(^) have the form 


ttQ 


arctan 


'^o(</) - — 

arctan[J/(20A)] ’ 

Uq = 2kl = const, A = 


[J/(20A)] ’ 
ttQ 


A > 0 , 
A < 0 , 

0 . 


n < cuo < 2Vt, 
2kl < ojq < 00 , 


System (11) is a system with two slow variables J,t and two fast variables ip,t. The 
existence and stability of stationary resonance solutions in such systems will be investigated. 
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4 Construction of averaged equations 


We will use the method of averaging on an inhnite interval (see [6]). The major part 
in all the problems related to the principle of averaging is in the changes of variables. 
The changes allow to eliminate fast variables from equations of motion within the given 
accuracy and thus separate slow motion from the fast one. 

Let Jpq be a solution of the equation 


where p, q are relatively prime integers. By making a change tl: = ip + we transform 
system (11) into 

^ = -4:euo{J)[f{t,T) + -fJuo{J)K^{ip + 

J =Uo{J) - - AeuoiJ)J~^[fit,T) +'yJuo{J)K^{p + Jz/f, J)]x (12) 

{-Jk{(p + ^ut,J))j. 

The point J = Jpq is resonant point in system (12). 

Assume that the resonance is non-degenerate, i.e., 

= UjQ^Jpq) 7^ 0. (13) 

We shall study the behavior of solutions of system (12) in the p = i/e-neighborhood of 
the resonance point Jpq. We make a change 


dujQ 

~dJ 



J — Jpq pZ 


and expand the right-hand side of system (12) in terms of the powers of /i. As a result, 
we obtain the system 


§ = pFo{t,T,ip,Jpq) + p‘^Fi{t,T,(p,Jpq)z + 0{p^), 

^ = fj,u’o{Jpq)z -h \p^uj'^{Jpq)z‘^ + /i^Go(t, r, (p, Jpq) + 0(/i^) 


where 


Fo(f,r,<p, Jpq) = -AuJo{Jpq)[f{t,T)+'yJpqUo{Jpq)K^{Lp+^Vt,J + P?)] X 
iZtpijp + pZ^^) Jpq)i 


Fi(t,r,(p, Jpq) = -Ajj |a;o(J)[/(t,r) + 7 Ja;o(J)«p(<p + fz/f, J)]x 
+ J)| 


.J=Jpq 


Cro{t, T, if, Jpq) Jpq) Jpq '^) T J Jpq^oi^Jpq) ^ipi^P T p^^) ^ 

( T p^^’ Jpq) Jpq^^jiSP T p^^) Jpq)' 


5 



System (14) contains only one fast variable t. We now make the standard change of the 
method of averaging in order to eliminate the fast variable in the right-hand side of system 
(14) up to the accuracy of the terms of order This change is sought in the form 

z = i + /iMi(?7, f, r) -h t, r)^, cp = t] + /i^n2(h, t, 't), 

where the functions Ui{ri,t,T), {i = 1,2), V 2 {ri,t,T) are periodic in t, r with period 27r/z/ 
and 27r/r accordingly. 

The change results in the system 


(15) 


ft = /^/o(h, + o{p^), 

f = t) + C>(/i^), 

where / o (? 7 , t ), / i (? 7 , t ), goirj^r) are dehned as the mean values over t: 

2nliy 

fo{v,T) = ^ f Fo{t,T,ri,Jpq)dt, 

0 

2n/v 

= I Fi{t,T,ri,Jpg)dt, 

0 

27r/iy 

go{v,r) = ^ f Go{t,T,g,Jpq)dt. 

0 

The functions Ui{g,t,T), {i = 1,2), V2{g,t,T) are dehned as periodic solutions in t with 
the zero mean value from the equations 


= Fi{t,T,g,Jpq) -uip{g,t,T)uJo{Jpq) -/i(^?,r). 


^ = Fo(t, r, r/, Jpq) - /o(? 7 , r), 

du2 _ 

~ 

= Go{t, r, g, Jpq) - ui{g, t, T)ujo{Jpq) - go{g, r). 

Functions /o(?7, t), fi{g, r), go{g, r) are smooth in g for considered perturbations. The 
function fo{g, r) is mean value over t of the function Fo(t, r, (p, Jpq). After averaging only 
a single term remains from the inhnite sum. Calculation of fo{g,T) will be demonstrated 
below (see (34)). Similar assertions hold for functions /i(? 7 ,r), go{g,T). 

System (15) at the time r is a singularly perturbed system in the following form 


= /o(h, + 0 {p‘^), 


dr 

= ^o{Jpq)i + \^^^o{Jpq))(^ + ^^9o{g, r) + 0(/i2), 


(16) 


5 Existence and stability of almost periodic solutions 

Let there exist periodic function go{T) such that 

/o(ho(T),r) = 0, 0 < 77o('r) < 27r. (17) 

In this case, the degenerate system derived from (16) at /i = 0 has the solution 

^ = 0, g = 77o(r). (18) 
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Linearizing the right-hand side of system (16) at /i = 0 on solution (18) yields the matrix 


Anir) = 


0 foviVo,r) 


uJq{Jj 


pqj 


0 


If 


(^o{Jpq)fov{Vo,r) > ao > 0, ao = const, r G (-cx),cx)), (19) 

then the matrix Ao(t) has real eigenvalues of different signs. In this case, as we know [6, 
chapter 8 ], in the system 

du ^ . , , 

/i— = Ao{t)u (20) 

for sufficiently small /i, the space of solutions U{fi) can be represented in the form 

u(p)^UM + U4^L). 

For the solutions M+(r,/i) G U+{fi) the inequality 

7+, 


|M+(r,/i)| < M+exp 


r - s 




|M+(s,/i)|, ( —CXD < S <T < OO) 


holds, and for U-{t, G the following inequality holds; 


|M_(r,/i)| < M_ exp 




|M+(s,/i)|, ( —CXD <T < S < Oo). 


Here M+,M_, 7 +, 7 _ are positive constants and | ■ | is a norm in It follows from an 
estimation of the solutions of system ( 20 ) that the solution of this system is unstable for 
sufficiently small fi if the space X_{fi) of the initial conditions of the solutions from U-{fi) 
is non-trivial. Hence, an inhomogeneous system 

d z 

L{fi) = fi— - Ao{t)z = fir), 

where /(r) is a periodic two-dimensional function, for sufficiently small fi has a unique 
periodic solution. This solution is represented as 


z{r,fi) = L {fi)f = - / K{T,s,fi)f{s)ds, 


fi 


where 


'y 

\K{t, s, fi)\ < Mexp ( — — \r — s\ ) (— cx) < t,s < oo), 

fi 


( 21 ) 


and M, 7 are positive constants. 

We transform system (16) using a change 


u = rf- Tfoir) 
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and write the obtained system in the vector form (z = 


dz 

jJ,— = Ao{t)z + H{z, t, T, /i). 
ar 

The components H{z, t, r, fi) have the form 

foiu + T]o, t) - fonivo, r)u + + ?7o, r)^ + 

+ figo{u + r]o, r) + 0{iZ). 

Evidently, the following inequality is valid 

<p{p), 


( 22 ) 


(23) 


where p(/i) ^ 0 as /x —)■ 0. The components of the vector function H{z, t, r, n) are 
differentiable at z in a sufficiently small neighborhood of the point (0,0). Therefore the 
following inequality holds 


-Z2I, 


(24) 


where pi(r, /i) —)■ 0 as r —)■ 0. 

The vector function II(z, t, et, fi) is an almost periodic function of variable t. The 
problem of almost periodic solutions of system (22) is equivalent to the problem of 
solvability of the system of integral equations 


z{t, = — / K{t,s, g)H{z,s, g)ds = Il{z, g). 
h" J 


(25) 


From inequalities (23), (24) follow that the successive approximations (see, for example, 

|ii|) 

«o(T,fi) = n(2,o), 2,(T,fi) = n(2j_i(r,fi),fi), j = i,2,... 

for sufficiently small fx converge uniformly on interval (—cx), cx)) to a unique almost periodic 
solution z*{T,fi) of the system of integral equations (25). The solution z*{T,g,) tends to 
(0,0) uniformly with respect to f as /i ^ 0. Hence, system (22), for sufficiently small fx, 
has a unique almost periodic solution z^{T,g). In its turn, system (16), for sufficiently 
small /i, has a unique almost periodic solution. Therefore, system (15), for sufficiently 
small /i, has a unique almost periodic solution. To investigate the stability of the almost 
periodic solution z^ir^jx) of system (22), we make a change 2 ; = z^,{t, fx) + y{T,fx) and 
obtain the system 


= Ai(p)y+ 


(26) 


where 


Hiiy, i), r, /i) = hf(^*(r, y) + |/(r, r, y) - H{z^{t, y),'i/j, r, y). 

The problem of the stability of the almost periodic solution z^:{T,y) is reduced to the 
problem of the stability of the zero solution of system (26). Taking into account the 



exponential estimates on the solutions of system (20), based on the Liapunov’s theorem 
of the instability by hrst approximation, we obtain that, for sufficiently small fi, the zero 
solution of system (16) is unstable. Hence, the almost periodic solution of system 

(22) is unstable. We shall state this result as a theorem applied to system (14). 

Let Jpq he a eonstant such that 

^oiJpq) ~ 0 ) ^oi'^pq) ^ 

Let there exists a periodic function rjoi^r) (0 < Pq^t) <2ti) such that 

fo{Vo{r),r) = 0 


and inequality (19) holds: 


(^oiJpq)foriiVo:r) > do > 0, do = coust, T G (-CX), oo). 

In this case, system (If), for sufficiently small n, has a unique unstable almost periodic 
solution. 

Hence, given the conditions of Theorem 1 are met in the p-neighborhood of the 
resonance point Jpq, there exists a unique unstable almost periodic solution. 

Now let, instead of inequality (19), the contrary inequality holds 

oJo{Jpq)for,{Vo,r) < ai <0, = const, -oo < r < oo. (27) 

If inequality (27) holds, the eigenvalues of the matrix Aq{t) for all r are purely imaginary. 
Now we need to consider averaged equations of higher approximations (we used only the 
hrst approximation in satisfying inequality (19)). 

We consider a narrower neighborhood of the resonance point Jpq. We make a change 
in system (15) 

f = puf) + lafoir) + r), 

V = ho(^) + 

where Pq^t) satishes equation (17), periodic functions Co('r) and Vo{t) are solutions of the 
equations 


= ^o('4g)^o(r) + 9o{Vo{r),T), 

= fop(Po(T),T)Vo(T) + fl(po(T),T)fo(r) + 

(Fi(t, T, Po{t), .Jpq)Ui{po{T),t, T))t + {Fop{t, T, r/o(r), Jpq)v2{po{T),t, T))t 

respectively. Here (•)* is mean value in t for hxed r. These equations can be solved using 
the inequalities (19) and (27). The functions u^{t,T) is dehned as periodic solution with 
zero mean value from the equation 

^ = Fop{t, T, po{t), Jpq)v2{po{T),t, r) + Fi{t, r, r/o(r), Jpq)ui{pQ{T),t, r)- 
{Fopit, T, Pair), Jpq)v 2 {po{T),t, T))t - {F^it, r, poir), Jpq)ui{poiT),t, r))*. 
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The change transforms the system into 


^ = na{T)v + + e(r)n^] + 0(/i^), 

^ = ^c{t)u + n‘^d{T)v + 0{n^), 


where 

a{r) = fov{Vo{r),T), 6(r) =/i(r/o(r), r), c(r) = |/o^^(r/o(r), r), 
d(r) = n;o( Jpg), e(r) = ^op(ho(T), r). 

Condition (27) in new notation takes the form 

a(r)c(r) < ai < 0. 


( 29 ) 


As was noted, it follows from this condition that the eigenvalues of the matrix of the 
hrst approximation are purely imaginary for all r. We reduce system (28) to "standard 
formi.e., to the form where the matrix of the hrst approximation is zero. A detailed 
description of reduction of system (28) is contained in [6, chapter 16]. 

We obtain the system 


dr 

dB 

dt 


i 6(r) + e(r) 
i 6(r) + e(r) 


h'jr) 

h(T) 

h(r) 


A + 

B + f2{A,B,x,T) +0{ii), 


(30) 


where 


h(r) = 


d{^) 
L «(^) 


1/2 


x(AAt) = - [ [-a(s)c(s)]^/^ds. 
0 


The functions fi{A, B, x, t), / 2 (A, B, x, t) are periodic in y with period 27i. They contain 
the terms not lower than quadratic in A and B. Analysis of the system (30) gives the 
following theorem for the system (14) (see similar Theorem 16.2 in [6]). (Below, the mean 
value of a periodic function /(r) is denoted by (/(t)).) 


Let a resonance point Jpq meet the conditions of Theorem 1. Let the inequality 

a{T)d{T) < (Ti < 0, t E {—oo, oo). 


holds. Finally, let the inequality 


(&(r) + e(r)) ^ 0 

holds. Then system (If ) in the e-neighborhood of the resonance point, for sufficiently small 
e, has a unique almost periodic solution. This solution is asymptotically stable if 


{6(r) + e(r)) < 0, 


and unstable if 


(6(r) + e(r)) > 0. 
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Theorems 1 and 2 apply to the study of resonant solutions in the vibro-impact system 
(11). The resonance point is dehned from equation 


Trhl 


q 

= -V. 


TT — arctan-^^ P 


(31) 


^0 ( Jpq ) 

From equation (31) follows that 

^o{.Jpq) > 0- 

To calculate remaining coefficients of the equation (30) we need to specify the function 
Let the function f{t,T) is dehned by the formula (8). 

To calculate a(r) = fopivo, t) we need to average Fo(f, r, ip, Jpq). The hrst term of this 
function is 

n 

(32) 


q 


- Auo{Jpq)E{T) sin(z/t + I3{t))k^{p + -ut). 


Since 


q 


K^{(p + Jj/t) = -UoiJpq) 

p TT 


k=l 


k sin k{ip + |z/t) 
-p 


, Hq r2[cao(t/pq] , 


we have to average summands in the form 


q 


E{t) sin(z/t + /9(r)) sin k{if + -z/f), /c = 1, 2,.... 

P 

It is easy to see that the mean value of function (32) will be non-zero if and only if 
g = 1, p = n (n = 1, 2,...). For g = 1, p = n it equals 


2E(r)z/2 
7rn{Q‘^ — z/2) 

The mean value of the second summand 


cos(np — /9(r)). 


q 


q 


-4'yJpqUJo{Jpq)K^{(p + -ut)K^{(p + -ut) 


equals 




pq 


1 + 


P P 


In calculating the mean value, we used the following equality 

sin“^ Trhlo = 1 + 

Consequently, the function po('r) is determined as the solution of the equation 


cos(np — /3{t)) = 


jJpqT^n 2 ,,2 


4i?(r)z/^ 


(fl^ - z/^) 1 + 


4fl2A2\ 


p 

'-'pq 


■) 


= Anir). 


(33) 
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Since An{T) —)■ cxd as n —)■ oo, we see that equation (33) can have solutions only for a hnite 
number of the values of n. If equation (33) has solutions for a given value n (|A„(r)| < 1), 
then these solutions are determined by the following formulas 

/3(r) arccosA„(r) 21ti 

= — ±-^-, / = 0, 

n n n 

Calculating the derivative of the function /o(? 7 ,t) at the points ?7oi('r) yields 

fov{VOl{T),T) = ~ (34) 

and, therefore, (34) has a positive sign at n points and a negative sign at n points. 
Straightforward calculations similar to the ones done earlier show that 

(6(r) +e(r)) = (/i(r/o,r) + go^{r]o,T)) < 0. 

Theorem 1 and 2 imply the following result. If the resonance point Ji„ is a solution 
of equation (31), then, if e is sufficiently small and |A„(r)| < 1, equation (11) has n 
solutions in the y/e - neighborhood of the resonance point, which are unstable resonant 
almost periodic in f, and n solutions, which are asymptotically stable resonant almost 
periodic in t in the e-neighborhood of the resonance point. 

The case when the force is determined by formula (9) is more simple. The functions 
/o(? 7 ,r), /i(? 7 ,t), go{v,T) do not depend on r. Thus ? 7 o(t) = rjQ = const. We obtain 

- 9) = (l + 

and 

9 arccosA^ 21% 

Voi = -± -^ +-, / = 0,...,n-l. 

n n n 

The variable r is contained only in terms of order 0{g?) in systems (15) and (28). We get 
a result similar to result obtained for the perturbation (8). 

6 Conclusion 

In this paper we have shown that in the considered vibro-impact system with one degree 
of freedom under small biharmonic perturbation with close and very different frequencies 
can occur stable and unstable almost periodic resonant modes. 

Our study is based on a combination of the method of averaging on inhnite interval 
[6] and the methods of singular perturbation theory [10,13]. The change of variables of 
method of averaging transforms the original system into a system with a smooth principal 
part. This system is then singularly perturbed in the slow time. Further study uses some 
ideas from the theory of singular perturbations. 

Similar methods can be used to study the resonant modes arising from the action 
of biharmonic perturbation in various problems of theory of oscillation described by 
equations with smooth and discontinuous coefficients. 
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